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Abstract. Ultra-cold atomic systems provide a new setting where to investigate the role of long- 
range interactions. In this paper we will review the basics features of those physical systems, in 
particular focusing on the case of Chromium atoms. On the experimental side, we report on the 

hi . observation of dipolar effects in the expansion dynamics of a Chromium Bose-Einstein condensate. 

(— I ' By using a Feshbach resonance, the scattering length characterising the contact interaction can be 

strongly reduced, thus increasing the relative effect of the dipole-dipole interaction. Such experi- 
ments make Chromium atoms the strongest candidates at present for the achievement of the strong 
dipolar regime. On the theoretical side, we investigate the behaviour of ultra-cold dipolar systems 
in the presence of a periodic potential. We discuss how to realise this situation experimentally and 
we characterise the system in terms of its quantum phases and metastable states, discussing in detail 
the differences with respect to the case of zero-range interactions. 
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L^; 1. INTRODUCTION: ULTRA-COLD ATOMIC SYSTEMS 

^ . The cooling and trapping techniques for neutral atoms developed in the last 20 years 

(N I have allowed the achievement of Bose-Einstein condensation [1], quantum degenerate 

^ ' Fermi gases [2], and the creation of several kinds of correlated systems. In all those 

^^ . cases, a crucial role is played by the statistics of the atoms, the interatomic interactions, 

the geometry of the system and their interplay. 

The first experiments were carried out with ultra-cold bosons in the weakly interacting 
P;" ■ regime, where the Gross-Pitaevskii equation provides a good description of the system. 

' Due to the very low temperature in those systems, interactions are usually s-wave and 

hence well described by a zero-range contact potential. Within the Gross-Pitaevskii 
^ ■ theory, a huge variety of phenomena has been described, like collective oscillations [3], 

interference [4] and coherence effects [5], non linear atom optics (4- wave mixing [6] 
and solitons [7]), superfluidity (sound propagation [8], scissor modes [9], and quantised 
vortices [10]). In polarised (single species) fermionic systems s-wave interactions are 
forbidden due to the Pauli principle, so that interactions are normally introduced by 
considering two different internal sublevels, mixtures of different species or p-wave 
interactions. 

By increasing the strength of the interactions (e.g. by enhancing the scattering length 
via Feshbach resonances [11]) or by changing the dimensionality of the system, corre- 
lations can be introduced. Striking examples in this direction were the observation of 
the superfluid to Mott insulator transition in optical lattices [12] and the Berezinskii- 
Kosterlitz-Thouless transition [13], the realisation of a Tonks-Girardeau gas [14], the 
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study of the BCS-BEC crossover in fermionic systems [15] and the creation of ultra- 
cold molecules [16]. 

A further possibility is to consider atoms interacting via long-range dipolar interac- 
tions [17]. This will be the topic of this paper, which summarises the talks given at the 
conference Dynamics and thermodynamics of systems with long range interaction: the- 
ory and experiments (Assisi, Italy - July 2007) by Chiara Menotti from ICFO, Barcelona 
and Tobias Koch from the group of Tilman Pfau at the University of Stuttgart. We will 
explain how long-range interactions arise in ultra-cold atomic systems, and in particular 
present the observation of weak and strong dipolar interactions in a degenerate sam- 
ple of Chromium atoms. Moreover, we will discuss the effect of dipolar interactions on 
atoms trapped in an optical lattice and point out the striking differences with respect to 
the case of zero-range interactions. The most recent experimental and theoretical results 
presented here are contained respectively in our papers [18, 19]. 

For the interested reader, extensive reviews on the physics of ultra-cold bosonic and 
fermionic ultra-cold gases can be found in [20, 21, 22, 23, 24]. 

2. INTERATOMIC INTERACTIONS 

Although ultra-cold gases have extremely low densities (typically below 10^^ cm~^), 
their properties are strongly influenced by interactions [20, 21, 22]. Usually, only the 
short range, isotropic contact interaction plays a role in quantum degenerate gases. 
However, another type of interaction, namely the anisotropic, long range interaction 
between dipolar particles, has attracted a lot of interest recently [17]. 

2.1. Contact potential 

In the ultra-low temperature regime relevant for quantum-degenerate gases, scattering 
occurs only in the s-wave regime, as the centrifugal barrier for other partial waves is 
much higher than the typical kinetic energies of the atoms. A consequence of this fact is 
that the real, complicated interatomic molecular potentials (which at long distances are 
essentially given by Van der Waals attraction —C(^/r^) can be replaced for most purposes 
by a simple, isotropic and short-range model potential proportional to the scattering 
length a of the atoms. This contact interaction reads 

AttPi a 
f/contact(r) = 5{v)^g5{v), (1) 

m 
where m is the atomic mass. Most of the interesting properties observed in quantum 
gases since the first observation of BEC in dilute atomic vapours in 1995 can be under- 
stood by using this simple interatomic potential. 

2.2. Dipole dipole interaction (DDI) potential 

Atoms or molecules having a permanent dipole moment (either magnetic or electric) 
interact not only via short-range potentials, but also via the dipole-dipole interaction. 
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Repulsive Attractive 

FIGURE 1. (a) Notations for the dipole-dipole interaction, (b) Dipoles placed side-to-side repel each 
other (c) Dipoles in a 'head-to-tail' configuration attract each other. 



The corresponding potential is 

J J I ^ Cddl-3cos20 

^dd(r) = — ^3 , (2) 

where Qd is the dipolar coupling constant (Cdd = /ioM^ for magnetic moments /i, 
Qd = ^^/£o for electric dipole moments J), and Q the angle between the direction 
joining the two dipoles and the dipole orientation (we assume here that all dipoles are 
aligned along the same direction z). The DDI is anisotropic (dipoles placed side-to- 
side repel each other, while dipoles in a head to tail configuration attract each other, 
see Fig.l) and long range (the 1/r-^ dependence implies for example that the scattering 
cross-section is not isotropic in the low-energy limit). 

To characterise the relative strength of the dipolar and contact interactions, it is 
convenient to introduce the dimensionless parameter 

Cdd = ■- .2 • C^) 

Miitra 

The numerical factors in £dd are chosen such that a homogeneous BEC with £dd > 1 is 
unstable (see Sect. 3. 2.1 below). For the alkalis usually used in BEC experiments, the 
value of Cdd is extremely small (for example, for ^^Rb, one has Cdd — 0.007), making the 
effects of the DDI negligible. 



2.3. Ultra-cold systems with long-range interaction 

There are several candidates to realize experimentally a dipolar quantum gas: 
molecules having a permanent electric dipole moment J, Rydberg atoms, which can 
have very large induced electric dipole moments, or ground state atoms having a large 
magnetic moment \x. In this section, we are going to briefly describe the main character- 
istics of those physical systems, and in the rest of this review, we will describe in detail 
the case of ultra-cold Chromium atoms, the only ones for which at present quantum 
degeneracy has been achieved. 
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2.3.1. Polar molecules 

Heteronuclear molecules in their ground state have a large electric dipole moment, on 
the order of one Debye (1 D ^ 3.3 x 10^^ C ■ m). Assuming that the order of magnitude 
for the scattering length is similar to that of atoms used in BEC experiments (typically 
around lOO^o, where uq is the Bohr radius), the corresponding value of £dd is on the 
order of 100, meaning that the properties of such a quantum gas would be dominated by 
the DDL 

To date, no quantum degenerate gas of polar molecules is available experimentally. 
Progress has been made recently in cooling of molecules (see ref. [25] for a review), but 
the densities and temperatures achieved so far are still orders of magnitude away from 
the quantum degenerate regime. A promising approach is to use Feshbach resonances 
in mixtures of ultra-cold fermions in order to create heteronuclear bosonic molecules. 
Those molecules created in a highly excited vibrational state must then be brought to the 
lowest vibrational state, e.g. by photoassociation. 



2.3.2. Rydberg atoms 

Another system in which large electric dipole moments can be achieved is given by 
Rydberg atoms. They are highly excited atoms, having large principal quantum number 
n and a dipole moment scaling like n^. The mutual interaction depends on the atomic 
states involved. By means of an electric field, one can mix states with different electron 
orbital angular momentum, such that the atoms acquire a dipole moment and can interact 
to first order via dipole interaction. Due to the very large dipole moments, this interaction 
is very strong and is felt over very long distances, of the order of many tens of microns. 

The state of the art in this field includes the observation of the blockade effect, 
which forbids more than one atom to be excited in a given region of space (defined 
by the distance where the interaction energy equals the linewidth of the excitation) 
and allows the production of an atomic ensemble with a single collective excitation. 
Moreover it includes the investigation of the resonant dipole-dipole interaction in an 
electric field (Forster resonance), collective behaviours and the real-time study of the 
dynamics of interacting pairs of Rydberg atoms, revealing the character of the long- 
range interactions. For more details, see e.g. [26] and references therein. 



2.3.3. Paramagnetic atoms with large magnetic moments ji 

Some atoms like Cr, Eu, Dy, have a large magnetic moment of several Bohr magnetons 
in their ground state. Only Cr has been Bose-condensed to date, and therefore the only 
quantum gas to display measurable dipolar effects is the Chromium BEC obtained in 
Stuttgart in 2004 [27] . Chromium has a magnetic dipole moment of 6/iB , and a scattering 
length of about IOOoq ("^o is the Bohr radius). This gives e^d — 0.16, which allows to 
observe a perturbative effect of the dipolar interaction. 
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3. BEC OF CHROMIUM 

In this section, we give a very brief overview of the experimental sequence used to obtain 
a ^^Cr BEC. More details can be found for example in [28]. 

The specific level structure of Chromium makes it possible (and necessary) to use 
novel laser cooling strategies to load atoms continuously into a magnetic trap [29]. 
After Doppler cooling in the magnetic trap, we get a cloud of 1.5 x 10^ atoms at a 
temperature of a few hundreds of /iK [30]. RF-induced evaporative cooling is then 
performed. However, dipolar relaxation from the low-field-seeking state I^Ss, m5 = +3) 
towards lower ms states prevents condensation [31]. 

The cloud (with 6 x 10^ atoms at about 20 /iK) is thus transferred into a crossed 
optical dipole trap (50 W at 1070 nm), and atoms are optically pumped to the high-field- 
seeking state pS3,m5 = —3), which is the absolute ground state; with a magnetic field 
of a few Gauss, dipolar relaxation is thus energetically suppressed. Forced evaporative 
cooling in the dipole trap then yields a pure condensate with up to 10^ atoms [32]. 



3.1. Feshbach resonances in Chromium 

Besides its large magnetic moment, Cr has another asset: 14 Feshbach resonances 
have been observed [33] for atoms in the state 1^83,^5 = —3). Close to such a resonance, 
the scattering length varies with the applied magnetic field B as 



a(5)=abg(^l-^-^J. (4) 

Here, Bq is the resonance position, A its width, and abg the background (non-resonant) 
scattering length. This yields the possibility of increasing e^d by making a approach 
zero. For this, one needs to control the magnetic field B with a precision much better 
than A. The broadest known Feshbach resonance of Cr lies at 5o = 589 G, and has a 
width of A = 1.4 G. The small value of A/Bq ~ 2.3 x 10^^ implies that one needs a 
good control of the field to tune a accurately (typically, a relative control of the field at 
the ~ 10^^ level is needed to control a at the ~ gq level). 



3.2. Dipolar expansion 

The most spectacular effect of the magnetic dipole-dipole interactions (MDDI) on 
the Cr BEC in an anisotropic harmonic trap appears in time of flight experiments. The 
aspect ratio of the cloud during expansion is modified by the MDDI and depends on 
the orientation of the atomic dipoles with respect to the trap axes. In this section, we 
first give a summary of the theoretical tools used to describe such experiments, and then 
describe our results. 
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3.2.1. GPEfor dipolar gases 

Pure contact interaction: a reminder. Weakly interacting BECs with pure contact 
interaction are well described by the Gross-Pitaevskii equation (GPE) for the order 
parameter i//(r, t) 

ih^ = -^^¥+ {V,n + g\¥\^) ¥■ (5) 

The non-linear term proportional to g accounts for the effect of interactions within the 
mean-field approximation. Note also that in the time-independent case, the left-hand 
side of the above equation has to be replaced by /ii//, with ji the chemical potential. The 
normalisation of i// chosen here is / | i//p = A^, where A'^ is the total atom number. 

A useful reformulation of the Gross-Pitaevskii equation is obtained by writing i// = 
V^exp(/5), with n the atomic density and S the phase of the order parameter, related to 
the superfluid velocity field by v = {h/m)VS. Substituting in (5) and separating real and 
imaginary parts, one gets the following set of hydrodynamic equations 

|^ + V-M=0, (6) 

the equation of continuity, and an Euler-like equation 

m^ + V — +,„ + V„--^)=0. (7) 



For the case of a uniform condensate (Vext = 0)^ one easily shows, by linearising the 
hydrodynamic equations around equilibrium, that the frequency o) and wavevector k of 
a harmonic perturbation are linked by the following dispersion relation, the Bogoliubov 
spectrum 



, jsn h^k^ 
m 4ot^ 

Dipolar interaction. To include dipolar effects, one just needs (as long as the DDI 
strength is not too high) to add an extra term to the mean-field potential g| V^P, namely 



^dd(r,0 = /lV^(r',Ol'f/dd(r-rVr'. 



(9) 



This extra term is thus non-local (due to the long-range character of the DDI) and makes 
it much more complicated to solve the GPE, even numerically (one faces now an integro- 
differential equation). For very large dipolar interaction (in the case, e.g., of dipolar 
molecules), this simple modification of the Gross Pitaevskii equation is not valid any 
more, as the separation of length scales on which it relies (namely, that the scattering 
length entering the contact interaction is determined by the intermolecular potential at 
relatively small distances, where the effects of the DDI are small compared to the van 
der Waals interaction) breaks down. In that case, the contact interaction depends on the 
strength of the dipolar coupling. 
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FIGURE 2. Density distribution for a non-dipolar BEC in an isotropic trap (left). The resulting dipolar 
potential Odd has a saddle-like shape, which tends to elongate the condensate along the magnetization 
direction (right). Figure taken from [27]. 



As a simple application, one can calculate the modifications to the Bogoliubov spec- 
trum induced by the DDI [34]. Using the fact that the Fourier transform of the DDI (2) 
takes the form 



t/dd(k)=Cdd(cos2a-l/3), 



(10) 



where a is the angle between k and the direction of the dipoles, and following the same 
method as in the previous paragraph, one easily shows that the excitation spectrum is 
now given by 



(O 




, + ^(3cos2a-i; 



+ 



4^' 



(11) 



and that, with the definition (3) for Edd, this implies that a dipolar uniform condensate is 
unstable for £dd > 1 • 



3.2.2. Thomas-Fermi solutions and scaling Ansatzfor the expansion 

Static Thomas-Fermi solutions. For pure contact repulsive interaction, when the 
atom number is increasing, the condensate size increases and the zero-point kinetic 
energy becomes smaller and smaller. The Thomas-Fermi approximation consists in 
neglecting the kinetic energy term in the time-independent GPE; this gives then a simple 
algebraic equation, showing that the density distribution has the shape of an inverted 
parabola. 

It is remarkable that this property remains valid if dipolar interaction is included. This 
is due to the fact that the dipolar mean-field potential ^dd(r) for a parabolic density 
distribution n(r) = | ^f{r) p is quadratic in the coordinates (having a saddle shape because 
of the anisotropy. Fig. 2) [35]. 

For the case of a spherically symmetric trap (and thus also a spherically symmetric 
density distribution for pure contact interaction), one can easily show that, to first order 
in Cdd, the effect of the MDDI is to elongate the condensate along the direction of 
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magnetization: it is energetically favourable to accommodate new particles close to the 
magnetization axis, where $dd(r) is minimum (see Fig. 2), thus causing an elongation 
of the condensate. It is possible to show that this behaviour is valid for anisotropic traps 
and for higher values of £dd. Note however that for a non-spherical density distribution, 
calculating the coefficients of the quadratic terms in 4>d(j is possible but beyond the scope 
of this paper. See [35] for details. 

Expansion: scaling Ansatz. For a pure contact interaction, and in the Thomas-Fermi 
approximation, namely Na/aho ^ 1 with aho = \/h/{inCo), there exists a very useful 
solution of the GPE [36]. It shows that the inverted parabola shape of the condensate 
is maintained upon expansion (after release from the trap), with a mere rescaling of 
its radii. The scaling parameters bi{t) (i = x,y,z) giving the radii Ri{t) = Ri{0)bi{t) are 
solutions of the ordinary differential equations 

^/- , ""ii^K {i^U,y,z}). (12) 

^i 11 ^■ 

je{x,y,z} 

Like for the static case, this can be extended to the case of dipolar interactions, since 
4>(id keeps a parabolic shape. The corresponding equations now read 

^' = u ""^K +n{bj};{cOj};eii) {i e {x,y,z}), (13) 

je{x,y,z} 

where / is a function of the scaling radii, trap frequencies, and dipolar parameter Cdd- It 
turns out that the elongation of the condensate along the magnetization direction remains 
valid during expansion. The reader is again referred to [35] for further details. 



3.3. Experiments 

3.3.1. MDDI as a small perturbation 

The first demonstration of an effect of the MDDI in a quantum gas came soon after 
the first realization of a Cr BEC, by measuring the aspect ratio of the BEC during time- 
of-flight for two different orientations of the dipoles with respect to the trap axes. The 
small value Cdd — 0.16 implied that the effect was only a small perturbation on top of 
the expansion driven by the contact interaction (see Fig. 3). 

3.3.2. Use of the Feshbach resonance: strong dipolar regime 

To go beyond this perturbative effect, we used the 589 G Feshbach resonance of Cr, 
in order to reduce a, and thus enhance £dd- We provide this field using the offset coils of 
the magnetic trap, with a current of about 400 amperes, actively stabilised at a level of 
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FIGURE 3. MDDI as a small perturbation in the expansion of a condensate. The aspect ratio is 
measured during expansion for two different orientations of the dipoles with respect to the trap axes. 
Figure taken from [27]. 



4 X 10^^ in relative value (peak to peak). The field is switched on during the evaporation 
sequence in the ODT, at a stage when the density is not too high, in order not to lose too 
many atoms by inelastic losses when crossing the Feshbach resonances. The rest of the 
experiment is performed in high field. After a BEC is obtained, we ramp the field close 
to the resonance in 10 ms, hold the field there for 2 ms, and take an absorption picture 
(still in high field) after 5 ms of time of flight. 
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FIGURE 4. (a) Measured scattering length a(B) across the Feshbach resonance at Bq = 589 G. The 
condensate images (i)-(iv) in the top of the figure clearly show a modification of the size (due to the 
reduction of a) and of the aspect ratio (due to enhanced dipolar effects) of the condensate, (b) Aspect 
ratio of the BEC after 5 ms of expansion, as a function of the measured dipolar parameter Edd- The solid 
line corresponds to the prediction of hydrodynamic equations in the Thomas-Fermi limit. Figure taken 
from [18]. 
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FIGURE 5. Aspect ratio of the BEC vs time of flight, (a) e^ = 0. 16, magnetic field along x, the aspect 
ratio is the same for the two trap configurations as expected, (b) e^^ — 0.16, magnetic field along z, one 
basically recovers the results of Fig. 3. (c) en = 0.5. d: e^ = 0.75, the inversion of ellipticity of the cloud 
is inhibited by the MDDl. Figure taken from [18]. 



From the density distribution, we measure the Thomas-Fermi radii of the BEC, and 
we infer the value of the scattering length (taking into account explicitly the MDDI 
interaction by solving Eqs. (13)). The measured a is shown on Fig. 4 (a). One can see 
clearly a five-fold reduction of a above resonance, corresponding to a maximal value 
of Edd — 0.8. On the sample absorption images of Fig. 4 (a), one clearly sees, when B 
approaches Bq + A, a strong reduction of the BEC size, due to the reduction of a, and 
thus of the mean field energy released upon expansion'. But one also clearly observes 
an elongation of the BEC along the magnetic field direction z. This change in the cloud 
aspect ratio would not happen for a pure contact interaction and is a direct signature of 
the MDDI. Fig. 4 (b) shows the aspect ratio of the cloud as a function of £dd, together 
with the theoretical prediction from (13). 

As an application of the tunability of Edd, we measured the aspect ratio of the BEC 
during expansion for two different orientations of the dipoles with respect to the trap 
axes. The effect of MDDI is now way beyond the perturbative regime (Fig. 5). For large 
enough £dd, one clearly sees that the usual inversion of ellipticity of the BEC during 
expansion is inhibited by the MDDI. 

Such tuning of the scattering length by means of a Feshbach resonance is a first 
step in the study of the properties of purely dipolar quantum gases (Cdd ^1). Future 
directions of study include e.g. the dependence of the stability of a dipolar BEC on the 
trap geometry, or the behaviour of dipolar quantum gases in optical lattices. 



For a pure contact interaction, the TF radius after expansion scales as (Na)^'^. 
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4. ULTRA-COLD ATOMS IN OPTICAL LATTICES 

The possibility of trapping and manipulating cold atomic gases with the degree of free- 
dom and precision described above allows to investigate a huge range of different physi- 
cal systems, where one has control on the system geometry, the interatomic interactions 
and the statistics of the atoms. One of the most fruitful fields of research both from the 
experimental and the theoretical point of view in the last years has been the study of 
ultra-cold atomic samples in optical lattices. 

Optical lattices are non dissipative periodic potential energy surfaces for the atoms 
created by the intersection of laser fields. The investigation of cold atoms in periodic 
potentials allows to reproduce problems traditional of condensed matter and solid state 
physics in a new setting, where a high degree of control is possible and where the 
Hamiltonian which governs the system is in general very close to some idealised one. 

One of the greatest achievements of the last years was the experimental observation 
of the superfluid to Mott insulator transition (for details, see Sect.4.1) [12, 37, 38]. 
In presence of long-range interactions, additional quantum phases are predicted [39, 
40, 41, 42], like the supersolid phase, presenting the coexistence of superfluidity and 
periodic spatial modulation, whose existence has not yet been unambiguously proven 
experimentally [43]. The search for this and other quantum phases makes cold atoms 
with long range interactions particularly appealing. 

As we will explain in the following, there are striking effects of long range interactions 
to be expected for ultra-cold atoms in optical lattices. Even if their existence strictly 
depends on the long-range character of the interactions, their observability might require 
a relative strength of the dipole-dipole interaction not too small compared to the zero- 
range one. For this reason the recent achievements with Chromium atoms presented 
in Sect. 3. 3. 2 are particularly important for the possible experimental realization of the 
systems that we are going to discuss in the following. 

First, we will introduce the main theory for ultra-cold atoms in optical lattices in the 
extensively studied and experimentally already realized case of point-like interaction. 
This will provide the background for the topic of our interest, namely the case of long- 
range interactions. 



4.1. Point-like interaction: superfluid to IVIott insulator transition 

As previously explained, in the most common cases, ultra-cold atoms interact via 
s-wave scattering, which can be in a very good approximation considered a point-like 
interaction. Then, the Hamiltonian of the system in second quantisation reads 

-PL j ^\r)^f{r) dr, (14) 
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where ^(r) is the field operator, Vextlr) the trapping potential, g the interaction strength, 
linked to the s-wave scattering length through g = Anh^a/m, and ji the chemical poten- 
tial, which fixes the total number of particles. 

In the presence of an optical lattice, the trapping potential is periodic Vextlr) = 
Y^n^n^^ ^i^^i'^^n/dn), where the index n runs over the dimensions of the lattice, and 
where dn is the lattice constant in the n-th direction. For lattices created by counterprop- 
agating laser beams of wavelength X , the lattice spacing is J = A /2. The intensity of the 
lattice potential Vn^^ depends on the intensity of the laser light. 

As it is well known [44], the spectrum of a single particle in a periodic potential is 
characterised by bands of allowed energies and energy gaps. For deep enough periodic 
potentials, the so-called tight binding regime is reached, where the first band takes 
the form E{q) = —2J'£„cos{qnd), being ^„ the quasi-momenta in the different lattice 
directions and J is the tunneling parameter between neighbouring wells. 

Provided the lattice is deep enough and for low enough interactions and temperature, 
the physics of the system can be well approximated by the one taking place in the first 
energy band. The counterpart of the energy eigenstates delocalised over the whole lattice 
(Bloch states) are the wavefunctions localised at the bottom of each lattice site (Wannier 
functions). In the tight binding regime is often convenient to use the Wannier description. 

In the single band approximation, the field operator can be replaced by its single-mode 
expansion 

x[r{r) = Y,Mr)d,, (15) 

with dj being the annihilation operator for one boson in the Wannier function w/(r) 
localised at the bottom of lattice site i. Plugging this expression in the second quantised 
Hamiltonian, one gets 

H = j;^Jw*{r)(^ + V,,,{r)\wj{r)drdJdj+ (16) 

+ 2 ^ w*{r)w*{r)wi{r)wm{r)dr dJd'^jdid,n-nY,^l^'' 

i,j,l,m i 

where in the chemical potential term we have exploited the orthonormality of the Wan- 
nier functions. 

For deep enough lattices, where the overlap beyond nearest neighbours can be ne- 
glected, and defining 

J = -|w;(r)(^|--Fyext(r))w,+i(r)Jr, (17) 

U = gJ\w,{r)\Ur, (18) 

one gets, a part from constant terms. 
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(ij) 



—ni[ni-l)-iini 



(19) 



with Hi = d-di being the density operator at site i. 

This is the famous Bose-Hubbard Hamikonian, extensively studied in condensed mat- 
ter physics. In optical lattices the Hamiltonian parameters can be accurately controlled 
by changing the strength of the optical lattice: ramping it up increases interactions due 
to a stronger localisation of the wavefunctions at the bottom of the lattice wells, and at 
the same time exponentially decreases tunneling. 

When tunneling is suppressed compared to interactions, this Hamiltonian presents 
a quantum phase transition between a superfluid phase, characterised by large number 
fluctuations at each lattice site, and a Mott insulating phase where each lattice well is 
occupied by precisely an integer number of atoms without any number fluctuations. The 
nature of this phase transition and the qualitative phase diagram can be inferred based 
on very simple arguments. 

At zero tunneling 7 = and commensurate filling (exactly an integer number n of 
atoms per well), the interaction energy is minimised by populating each lattice well by 
exactly n atoms. Energy considerations can tell which is the range of chemical potential 
/i at which the filling factor n is the most energetically convenient: 

E{n) = —n{n — 1) —/in, 

E{n)<E{n+\) and E{n) <E{n-l)^ 

{n-\)U <il<nU. (20) 

This state with precisely integer occupation at the lattice sites is called Mott insulating 
state. Since a particle-hole excitation at 7 = costs an energy AE = U equal to the 
interaction energy, the Mott state is the lowest energy state at commensurate filling (see 
Fig. 6). For a tunneling 7 different from zero the energy cost to create an excitation 
decreases thanks to the kinetic energy favouring particle hopping. However, for large 
interactions and small tunneling, the gain in kinetic energy (~ 7) is not yet sufficient 
to overcome the cost in interaction energy (~ U), which leads to the existence of Mott 
insulating states also at finite tunneling. For large enough tunneling, instead, particle 
hopping becomes energetically favourable and the system becomes superfluid. The 
regions in the 7 vs. /i phase diagram where the Mott insulating state is the ground state 
are called Mott lobes. For non commensurate filling, there are extra atoms free to hop to 
site to site at no energy cost, so that the phase of the system is always superfluid. The 
superfluid phase at non commensurate densities survives down to 7 = for ji/U = [p] 
where the symbol [p] indicates the integer part of the density. 

Due to the finite energy cost required to add or remove one particle, the Mott phase 
is gapped and incompressible, while in the superfluid regions the gap vanishes and the 
system is compressible. 

In Eq.(20) we have identified the boundaries of the Mott lobes at 7 = 0. In order to 
find the shape of the lobes at finite 7 more sophisticated calculations are required. In 
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FIGURE 6. Sketch of a particle-hole excitation in the « = 1 Mott phase. The energy cost for having 
two particles in the same site is U and the energy gained by hopping is ^ /. Their interplay determines 
whether the ground state is insulating or superfluid. 



particular there is no exact analytical method which allows to calculate them. A part 
from the mean-field approximation [37, 45, 46], which is discussed later, high order 
perturbative strong coupling expansions can be performed [47]. Exact numerical results 
can be obtained using Quantum Monte Carlo techniques [48]. 

In the next sections, we introduce the mean-field approximation and the perturbative 
method that we use to calculate the lobes. It works only qualitatively in one dimension 
and work better and better in larger dimensions. This will provide further insight on the 
conditions required for the insulating lobes to exist and proves very useful to understand 
the analogies and differences in the case of long-range interactions. 

In the non uniform geometries available in the experiments, where a harmonic poten- 
tial is superimposed to the optical lattice, the trapping potential is taken into account in 
the Bose-Hubbard Hamiltonian through a site-dependent chemical potential iii = ji— Vt, 
being Vj the harmonic potential at site i. The system is characterised by alternating shells 
of Mott and superfluid phases [38, 49, 50]. In this paper we will concentrate on the uni- 
form system and not discuss the trapped system in more detail. 



4.1.1. Mean-field decoupling and Gutzwiller Ansatz 

Starting from the Bose-Hubbard Hamiltonian (19), one can perform the so-called 
mean-field decoupling [37, 45, 46]. We define the order parameter (a) = (p and proceed 
with the following substitution 



~ + ~t * 



obtaining for the hopping term 



{ij) L 



^t[^J-^)+^jU-\]+^]^J 



(21) 



(22) 



Based on this result, the full Hamiltonian can be rewritten with no approximations as 



(ij) 



<?/■ I ^j - ^) + <?'■ f ^J - % ) + ^l^j 
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where the term a]dj is supposed to be small. 

The full Hamiltonian in the mean-field approximation is just given by the sum of the 
mean-field Hamiltonians on different sites H'^^ = Y^iHf^, where 



H 



MF 



-J 



^i\a 



f)... 



2 



+ —ni{ni - I) - lini 



(24) 



and (pi = Y.{j}i 9j takes the neighbouring wells into account at the mean-field level. 

In the homogeneous system, all order parameters (pi = (p are equal. For each set 
of parameters t/,7,/i, one can find the value of (p which minimises the ground state 
eingevalue (at T = 0). A vanishing order parameter is the indication of a Mott insulating 
phase, while an order parameter different from zero, which can be only produced by 
number fluctuations at the lattice sites, is the signature of a superfluid state. The mean- 
field solution found in this way is completely equivalent to the one provided by the 
famous Gutzwiller Ansatz. This corresponds in writing the state of the system as a 
product over the different lattice sites of single-site wavefunctions 



M^)) = ULfn\^)M^ 



(25) 



where / is the site label and n indicates the number state. The coefficients /„ are called 
Gutzwiller coefficients. 

The solution obtained by the Gutzwiller Ansatz is intrinsically mean-field. On the 
other hand it allows to extract many important information and study the dynamics of the 
system. In fact, through a time dependent variational approach based on the Lagrangian 



(4) 1 4)) - (4) 1 4)) 
one obtains the dynamical equations for the Gutzwiller coefficients 



(26) 



.df;i 



{') 



dt 



-J 



(PiV^if:\+(ptV^^UnU 



+ 



—ni[ni-l)-lini 



f!f^. (27) 



By solving those equations in imaginary (t = it) time, one accesses the stationary states 
of the system, while in real time the dynamics of the system can be investigated also in 
the case of time-dependent Hamiltonian parameters [51]. 



4.1.2. Mean-field perturbative approach 

The mean-field perturbative approach allows to find the boundaries of the lobes in 
an analytical way. The main approximation underlying this method is the mean-field 
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approximation described in the previous section. For the purposes of this section, we 
separate the Hamiltonian in two parts 



Hj 



—ni[ni-\)-lini 



-^L 



'pri^/-f)+<pJ^J-Y 



(28) 
(29) 



where we remind that 9, = (a,) is the order parameter and (pi = £/y\ . (pj is the sum of the 
order parameters at sites neighbouring to a given site i. 

The term Hj replaces the non-local tunneling term present in the full Hamiltonian. 
Assuming that the order parameters % are small, we rewrite Hj at first order in the order 
parameter as 



Hj ^ -7£ 



(pi Ui + (piul 



(30) 



We use this approximate expression in the calculation of the partition function 



^ 



Tr 



-pHo 



j8 
- / Tr 

'0 



-{p-^n 



Hje 



-THo 



dT + 



Tr 



,-m 



(31) 



where we have used that Hj is first order in the creation and destruction operators and 
neglected second order terms. In the limit /3 — ;> 0°, corresponding to zero temperature, 
the partition function becomes 



3r=^o = e 



-pEo 



(32) 



being Eq the energy of the ground state. In the calculation of the order parameter instead, 
the first order term in the creation and destruction operators also contributes, in the form 



(Pi 



Ma 



Tr 



a,e-(^-^)^"//ye-^«» 



= y(p/e^^" / Tr 
Jo 






dx = 



(33) 



where we have used the definition of/// at first order in the order parameter, the fact that 
the expectation value of two destruction operators over eigenstates of //q vanishes, and 
the fact that non local averages vanish in the mean-field approximation. For the trace, we 
need to take into account all the matrix elements of the kind {^\die''^^~'^^^^d\e~'^^'^\^) . 
The two creation and annihilation operators in this expectation value produce 
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n;- + 1 ) exp 



-m 



l*+i,0 



exp 



TiHJ^^-H, 



■l*+i,') 



(34) 



where |4>+ 1, /) is the notation we use for the state which is obtained from state |4>) by 

creating one more particle at site i and n- is the number of atoms at site i in state |$). 
The integration over T in the range T G [0, /3] , in the limit /3 ^ oo and the multiplication 
by e^^" give 



.;*> + ! 



exp 



PiEo-H, 



*+i,i-) 



exp 



PiEo-H. 



*) 



H. 



*) 



H, 



*+!,/) 



H, 



*) 



//, 



l*+i,0 



(35) 



■I*) 



Since by definition of ground state Eq — Hq < for all states |$), this result converges 
to a non zero result if and only if either 



1) 1$) is the ground state or 

2) |4> + 1 , z) is the ground state, i.e. |4>) is the state 

obtained by removing one particle from the ground state at site i. 



(36) 



Namely, only two terms contribute to the trace (33): the one arising from the ground state 
{GS\die~^^l^~''''i^°dj e~'^^^\GS) and the one arising from the states which are obtained 

from the ground state by removing one particle {GS — 1 , /|a,e^(^^^^^"aje^^'^f' | GS — 1 , z) . 
Hence, for a Mott state with integer well occupation equal to n,, the order parameter 
reads 



(Pi = 9iJ 



m + l 



+ 



Uui-ll jl~U{ni-l) 



(37) 



In a uniform system in the presence of point-like interaction (p, is uniform and the 
solution can be found analytically. Given that 9 = z^, being z the number of nearest 
neighbours, the order parameter (p = Q unless 



\-zJ 



n + \ 



+ 



Un-jl jl-U{n-\] 



0, i.e. zJ — 



{Un-^){^-U{n-\)) 



(38) 



This equation has solutions (for zJ > 0) only for chemical potentials in the range 
n — l< ji/U < n, as already anticipated in Eq.(20). 

Expression (38) of /i as a function of 7 defines the boundary of the Mott lobe at filling 
factor n. Outside of the lobes the order parameter becomes different from zero, indicating 
the arising of the superfluid phase. Then, this treatment, valid at 1st order in (p, is not 
valid anymore to predict the value of the order parameter. However, in the framework of 
the mean-field approximation, the expression for the lobes boundary is exact. Moreover, 
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most important, the prediction about the existence of the lobes, i.e. the lobes boundaries 
at 7 = 0, are exact also beyond the mean-field approximation and confirmed by exact 
numerical methods. 

For the comparison that we will make later with the case of long-range interaction it is 
useful to state the condition for the existence of the lobes as follows: given an insulating 
state (classical distribution of atoms in the lattice), this is stable in a given range of 
chemical potential at 7 = if the energy increases by adding or removing a particle at 
any site of the lattice, as depicted in Fig. 7. 







FIGURE 7. Sketch of a Mott insulating n = 1 phase with energy E{n ~ I) (left); configuration where 1 
atom has been removed, with energy E (n = I , — I , i) = E {n = I) + ji (center); configuration where 1 atom 
has been added, with energy E{n = l,+l, i) = E{n = I) — ^ + U (right). For < ^ < C/, the Mott « = 1 
phase is the ground state. 



4.1.3. Phase diagram 

Starting from Eq.(38) we can built the phase diagram of the system in the mean-field 
approximation, as shown in Fig. 8. Inside the lobes the system in in the Mott insulating 
phase. For a given tunneling parameter J and chemical potential /i the energy for adding 
or removing a particle from the system is respectively given by the distance to the upper 
and lower lobe boundaries (at constant /). The sum of those energies (width of the lobe 
at constant J) gives instead the energy for a particle-hole excitation conserving the total 
number of atoms. 

At the lobe boundary a quantum phase transition to the superfluid phase takes place. 
The tip of the lobe is a special point, because there the phase transition happens at 
commensurate filling factor n and is purely due to the effect of increased tunneling. At 
constant non commensurate density the system is always superfluid. For 7—^0, the lines 
of constant density accumulate between the Mott lobes, at integer values of ji/U = n 
where the energies for having n or n + l particle per site are equal. 

5. LONG-RANGE INTERACTIONS IN OPTICAL LATTICES 

We now consider dipolar atoms in a 2D optical lattice. In present experiments, usually 
2D geometries are created as a series of pancake traps by means of a very strong 
ID optical lattice in the perpendicular direction, which provides strong confinement 
and completely suppresses tunneling in one direction. In the presence of long-range 
interaction, in order to consider each layer to be isolated, one should not only suppress 
tunneling but also reduce the interaction between the different layers. For our purposes, 
where each layer contains a 2D lattice, one should make the distance between the 
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)i/U 




FIGURE 8. Phase diagram in the mean-field approximation. The boundaries of the Mott lobes are given 
by Eq.(38). The upwards and downwards arrows describe respectively the energy of a particle and hole 
excitation, as described in the text. Their sum is the energy of a particle-hole excitation. 



different layers larger than the lattice spacing in the 2D plane. This can be achieved 
by creating a 2D lattice with two pairs of counterpropagating laser beams and the extra 
ID lattice in the perpendicular direction by two laser beams intersecting at a given angle 
6, as depicted in Fig.9, in order to increase the lattice spacing in the third direction to 

dw^{?i/2)/sini9/2). 
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FIGURE 9. Two-dimensional lattice potential obtained by the intersection of four counter-propagating 
laser beams. For a wavelength X the lattice spacing is given by d^o ~ X/2 (left); extra one-dimensional 
lattice in the perpendicular direction, obtained by the intersection of two laser beams at and angle 9; the 
lattice spacing is given by dio = {X/2)/sin{6 /2) (right). 
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We consider an orientation of the dipoles perpendicular to the planes of the 2D lattice. 
This implies a repulsive interaction in the plane and an attractive interaction between 
different layers. The effect of this intra-layer attraction has been studied in [52]. We 
expect that the attractive character of the interaction in the perpendicular direction might 
even help in reproducing in the different layers the distribution of atoms that we find to 
exist in the single plane. However, to draw exact conclusions one should devote to this 
topic further studies. In the following we are going to restrict ourselves to the study of a 
single lattice plane. 



5.1. Extended Bose-Hubbard model 

In the presence of long-range interaction, an extra term has to be added to the Hamil- 
tonian of the system, which in second quantisation reads 



H = |i?rt(r)(^|- + yext(r)ji?r(r)Jr+||i^t(r),^t(r),^(r),^(r)Jr + 

J Jxir'^{r')xif\r')Uid{r-r')xif{r)xif{r)dvdr'-^Jxif'^{r)Yir)dr, (39) 



+ 



where t/dd(r) is the dipole-dipole potential defined in Eq.(2). 

In the single band approximation, the generalised Bose-Hubbard Hamiltonian in the 
presence of long-range interaction becomes 



H 



(ij) 



—ni[ni-l)-iini 






(40) 



where £ is the distance connecting the two optical lattice sites i and j. The sum over 
the distance £ is cut-off at a certain nearest neighbour. In our calculations we usually 
considered up to the Ath nearest neighbour, as shown in Fig. 10. Longer interaction 
ranges have a crucial effect on the regions of the phase diagram describing low density 
of particles and holes [53], but we are not going to discuss this here. 
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FIGURE 10. Representation of the first four nearest neighbours in a 2D optical lattice. 
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The on-site interaction parameter U is now given by two contributions: one is, as 
before, arising from the i'-wave scattering U^ = ATzfP'a/m j n^{r)d^r, and the second one 
is due to the on-site dipole-dipole interaction C/dip = ^/(^Tl) jUdd{q)n^{(i)d^(i, being 
Udd{(l) and n{q) the Fourier transform of the dipole potential and density, respectively 
[54]. 

Because of the localisation of the wavefunctions at the bottom of the optical lat- 
tice wells, the long range part of the dipole-dipole interaction U^ is in a very good 

approximation given by the dipole-dipole interaction potential at distance i, Uj — 
(Cdd/4;r)[l — 3cos^(07')]/£-^ multiplied by the densities ni and rij in the two sites, where 
the quantity Cdd, proportional to the dipole moment squared, has been defined Sect. 2. 2 
and 0-^ is the angle between I and the orientation of the dipoles (see Fig. 1). 

The ratio between the total on-site interaction U = Us + f/dip and the nearest neighbour 
dipolar interaction Uf^j^ determines much of the physics of the system. As we mentioned 
previously, in order to make the effects of long-range interactions observable it might be 
necessary to have a not too small ratio U /Un^. It can be varied by tuning the strength 
and the sign of the on-site dipole-dipole interaction t/dip by changing the anisotropy 
of the Wannier functions at the bottom of the lattice sites, as depicted in Fig.ll. In 
standard experiments with ^^Cr, U /Umn ~ 400 (for £dd = lloi^ilBfrn/ {llnh^a) ^ 0.158 
and spherical localisation at the bottom of the potential well at 5 = 20£'^, where Er is 
the recoil energy at A = 500 nm). Using a Feshbach resonance to change the 5- wave 
scattering length, as recently demonstrated with Chromium atoms [18], U /Umn can be 
virtually tuned down to zero. 






Ump<0 U,,^=0 U,,^>0 

FIGURE 11. On-site dipole-dipole interaction J/dip depending on the anisotropy of the wavefunction 
at the bottom of the lattice wells. For vertically pointing dipoles, the dipole-dipole interaction is mainly 
attractive for cigar-shaped wells (left), vanishes for spherical wells (center) and is mainly repulsive for 
pancake-shaped wells (right). 

Analogously as before, a time dependent variational principle leads to the dynamical 
equations for the Gutzwiller coefficients which now include a contribution due to the 
long-range part of the interaction 
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The imaginary time evolution, which mimics dissipation in the system, converges un- 
ambiguously to the ground state of the system for the Bose-Hubbard model in presence 
of on-site interaction only. In the presence of long-range interaction it shows a strikingly 
different behaviour and converges often to different configurations, depending on the ex- 
act initial conditions. In this way, we clearly get a feeling of the existence of metastable 
states in the system. In the real time evolution, their stability is confirmed by typical 
small oscillations around a local minimum of the energy. 



5.7.7. Conditions to have an insulating lobe 

The existence of metastable stationary states can be confirmed using the mean-field 
perturbative approach introduced in Sect.4.1.2. In this section we will generalise it 
to include long-range interaction and we show that, beyond the ground state lobe, 
insulating lobes for the metastable states exist [19]. 

The first important difference is that in the limit /3 ^ oo, the system can populate 
any local minimum of the energy and not necessarily the ground state, as it was before. 
Hence, the partition function becomes 



ay — ay 

ZC = ZCfL 



MS 



,-i6£* 



(42) 



where MS indicates any of the metastable states, whose existence we are for the moment 
assuming and going to demonstrate in the following. 

The expectation value of the order parameter reads, similarly to before 



VTS-^/^^* 



(Pi = JipicP^'^' / Tr 
Jo 



P 



^.g-(/3-T)//o^tg-T//o 



dx, 



(43) 



where Hq now includes also the dipole-dipole interaction 



^0 = 1 



U 



—ni{ni-l)-^ni + Y,U^ 



Hi [rii- 



im) 



(44) 



On the other hand, the tunneling part Hj has not changed with respect to the case of 
zero-range interaction. 
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It is important to stress at that point that in order to investigate the stability of a given 
metastable state, the trace is performed on the subspace of states differing from it only by 
small perturbations. Specifically this means that we take into account only those states 
which differ from the metastable state by adding or removing one particle. 

The derivation follows the one presented in Sect.4.1.2. The important analogies and 
differences are summarised here: 

a) the ground state is replaced by any metastable state and only small perturbations 
around that metastable state are considered; 

b) the assumption of metastability of a given configuration has to ensured by checking 
that all the states obtained from the metastable state \MS) by adding or removing one 
particle at any of the lattice sites j have higher energy than Ems, i-e. 



Ems - ^r < 0' for |4>) = |M5 ± 1 , j) Vj; (45) 

c) the only non vanishing contributions arise, analogously to Eq.(36), from the terms 
where 



1) 1 4>) is the metastable state or (46) 

2) |<I> + 1, /) is the metastable state, i.e. |4>) is the state 

obtained by removing one particle from the metastable state at site /. 

These are exactly the same criteria as for on-site interaction only, a crucial difference 
being that now the atomic distribution in the lattice sites might be not uniform so that 
the effect of adding and removing one particle at all lattice sites has to be considered to 
ensure condition (45) above. Taking the case of dipole-dipole interaction explicitly into 
account, one gets that for all i the conditions 

Uin,-\)<^-V^;i^<Uni (47) 

have to be satisfied, implying for the boundary of the lobes at 7 = the values 

M„„-„ = max [f/ (n,- - 1 ) + VJ.^] , (48) 

llma^ = mm Uui + V^Jp , (49) 

where VJ is the dipole-dipole interaction of 1 atoms at site i with the rest of the lattice. 

Clearly in the absence of dipole-dipole interaction, these conditions reduce to the ones 

written in Eq. (20). 

Those stability conditions immediately lead us to draw two important conclusions: 
i) configurations with no integer filling factor can be stable as soon as one nearest 

neighbour is included in the interaction (see Fig. 12(a)); 
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FIGURE 12. (a) Sketch of a checkerboard Mott insulating phase with energy E{GS) (left); config- 
uration where 1 atom has been removed with energy E{GS, — l,i) = E(n ^ I) + ji (center); configu- 
ration where 1 atom has been added with energy E{n = 1,+1,/) = E{n = 1) — /X +AUnn (right). For 
< ji < A-Unn, the Mott « = 1 phase is the ground state, (b) Sketch of a metastable Mott insulating 
phase at filling factor 1/2 with energy E{MS) (left); configuration where 1 atom has been removed with 
energy E{MS, — 1, /)=£'(« = 1) +/i — Unn (center); configuration where 1 atom has been added with 
energy E{n = 1,+!,/) = E{n = 1) — /i +3Unn (right). For Unn < l-i < 3Unn, this metastable state is 
stable. For the sake of simplicity, in those examples we have considered only one nearest neighbour in the 
dipole-dipole interaction. 



ii) for a given chemical potential there exist configurations at higher energy compared 
to the ground state which also fulfill the stability condition (see Fig. 12(b)). We call such 
configurations metastable states. 

The explicit expression for the order parameter now reads 



<P/ 



9,7 



rii + l 



rii 



Um-^ + v'f 



dip 



[/(n,-i; 



M+Ki;, 



(50) 



Equation (50) is a system of linear equations in the variables (p, described by the matrix 
M: M(p = 0. When det{M) ^ 0, this equation only allows the trivial solution q>i = V/, 
implying that the system is in a Mott insulating phase. On the other hand, if det{M) = 0, 
it is possible to have (pi ^ 0, implying that the system has become superfluid. 

Due to the non homogeneity of the solution in the presence of long-range interaction 
this is a system of linear equations whose dimension depends on the size of the lattice. 
In practise we can access with this method lattice sizes of the order of 20-30 lattice sites 
per direction. 



5.1.2. Phase diagram for long-range interactions 

The phase diagram of the system is shown in Fig. 13. As predicted by the perturbative 
mean-field approach presented in the previous section, in the presence of long-range in- 
teractions, there are insulating lobes corresponding both to ground and metastable states 
and to integer and non integer filling factors. In particular non uniform distributions of 
the atoms are allowed. 
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FIGURE 13. Phase diagram for weak and strong dipole-dipole interaction and interaction range up 
to the 4th nearest neighbour: U /Unn = 20 (a) and U/Unn = 2 (b). The thick lines are the ground state 
lobes, found (for increasing chemicals potential) for filling factors equal to all multiples of 1/8. The thin 
lines are the metastable states, found at filling factors equal to multiples of 1/16. Some of the metastable 
configurations at filling factor 1/2 (I to III) and corresponding ground state (IV). Empty sites are light and 
sites occupied with 1 atom are dark. Figure from [19]. 



The phase diagram in Fig. 13 has been calculated for an interaction range of four 
nearest neighbours, as shown in Fig. 10. For this specific range of interaction, the lowest 
filling factor allowed is 1/8. 

By increasing the relative strength the long-range interaction with respect to the on- 
site one, the lobes corresponding to fractional filling factors become comparable to the 
ones relative of the standard Mott phase with commensurate filling (cfr. Fig. 13(b) with 
Fig. 13(a), where the upper lobe for n = 1 is only partially shown). 

This phase diagram is confirmed by the imaginary and real time evolution of the 
system: depending on the initial conditions the imaginary time evolution can converge to 
the metastable configurations, while in the real time evolution, their stability is reflected 
into typical small oscillations around a local minimum of the energy. In the next section 
we discuss the lifetime of the metastable states, connected to the possibility of tunneling 
between different local minima of the energy landscape. 
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5.2. Stability of the metastable states 

We study the stability of the metastable states with a path integral formulation in 
imaginary time [55], which can describe the tunneling below a potential barrier (instan- 
ton effect). The path integral representation for the propagator in imaginary time takes 
the form 



N 



^[^\^]Y[{^j 



\e-"^'\^i 



t)^ /'^[$*,4>]e--^^^(^'*'*)^'^ (51) 



which gives the following expression for the Lagrangian Jf = (^|4>) —H. With the 
substitution 4> = (x + ip)/^/2 and 4>* = (x — ip)/^/2, the Lagrangian becomes ^ = 
—ipx + H. This can be put in the canonical form ^ = PX —H{X,P) by defining the new 
coordinates X = x and P = dJ^/dX = —ip. In those new coordinates, the Hamiltonian 
governing the dynamics in imaginary time is 



H{X,P) = -H{p{P),X) = ^ - V(^)- 



(52) 



For this reason, the dynamics in imaginary time, describing the instanton effect between 
two potential wells, is said to happen in the inverted potential —V{X), as depicted in 
Fig.l4. 





FIGURE 14. Tunneling of a particle in a double-well potential (left); in imaginary time, this tunneling 
event (instanton) is described by the dynamics in the inverted potential (right). 

We exploit this formalisms to describe the tunneling between the different metastable 
configurations found in the previous section. In order to reduce our problem to the 
dynamics of a single variable X and its conjugate momentum P, we make a variational 
Ansatz on the Gutzwiller wavefunction which parametrises the population of the lattice 
well from 1 to and viceversa. 



FIGURE 15. Transition from a checkerboard configuration to its opposite one, where full sites are 
replaced by empty sites and viceversa. 
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Here we present in particular the results for the transition between one configuration 
(in the specific case, the checkerboard) and its opposite one, where full and empty 
wells have been interchanged, as shown in Fig. 15. For this transition, the tunneling time 
is cOqT ^ exp[5o], where ft\) is the typical oscillation frequency around the minimum 
of the energy. The tunneling time diverges for 7^0 and roughly scales like cOqT ^ 
exp[Nsh] exp[—NshJ/J] for J/J> 0.3, as shown in Fig. 16. 




FIGURE 16. Action ^o for the transition described in Fig. 15, calculated with a variational Ansatz for the 
path integral formalisms in imaginary time. The lifetime of the tunneling event is given by COq T w exp[5o] . 
Ns is the number of sites which invert their population in the transition and / is of the order of the tip of 
the insulating lobe. 

This analysis suggests that the metastable configurations are very stable when many 
sites must invert their population to reach another metastable state. However, one should 
not forget that especially in larger lattices two metastable configurations might differ just 
by the occupation of few lattice sites. This, and the corresponding small energy differ- 
ence, should be carefully taken into account in a realistic analysis at finite temperature. 



5.3. Initialisation and detection 

Very important issues are the initialisation and detection of the atomic states in the 
lattice. One can use superlattices in order to prepare the atoms in configurations of 
preferential symmetry. This idea is presently pursued by several experimental groups 
[56]. We have checked that the presence of defects is strongly reduced when a local 
potential energy following desired patterns is added to the optical lattice. Note that the 
configurations obtained in such a way will also remain stable once the superlattice is 
removed, thanks to dipole-dipole interaction. 

The spatially modulated structures created in such a way can be detected via the 
measurement of the noise correlations of the expansion pictures [42, 57, 58] : the ordered 
structures in the lattice give rise to different patterns in the spatial noise correlation 
function 



C{d) 
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equal to the modulus square of the Fourier transform of the density distribution in the 
lattice {ujoF is the density distribution after time of flight, while n.]^ is the density 
distribution in the lattice). Such a measurement is in principle able to recognise the 
defects in the density distribution, which could be exactly reconstructed starting from the 
patterns in the spatial noise correlation function. For the moment, the signal to noise ratio 
required for single defect recognition is beyond the present experimental possibilities. 
However, averaging over a finite number of different experimental runs producing the 
same spatial distribution of atoms in the lattice, a good signal can be obtained. In Fig. 17, 
we show the noise correlations for the metastable configurations at filling factor 1/2 
shown in Fig. 13. 
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FIGURE 17. Spatial noise correlation patterns for configurations (I) to (III) in Fig. 13, assuming a 
localised gaussian density distribution at each lattice site. Figure from [19]. 

Presently we are studying the possibility of transferring in a controlled way those 
systems from one configuration to another [53]. This, together with the capability of 
initialising and reading out the state of the lattice, might make those systems useful for 
applications as quantum memories. 

6. CONCLUSIONS 

In this paper, we have given a review of some aspects of the physics of ultra-cold atomic 
gases interacting via a long-range dipolar potential. On the experimental side, we have 
presented an overview of the state of the art of the experiments, starting from the first ob- 
servation of dipolar effects in a Chromium Bose-Einstein condensate to the most recent 
experiments demonstrating strong dipolar interactions in those systems. The description 
of those experimental achievements has been accompanied by the explanation of the un- 
derlying theoretical models. On the theoretical side, we have presented our recent results 
on dipolar gases in optical lattices, describing the theoretical framework and explaining 
in a detailed way the difference brought in by long-range interactions compared to the 
case of zero-range interactions. In particular, we have pointed out that long-range in- 
teractions introduce a huge number of metastable states in the system, which are not 
present in the case of zero-range interaction. We have motivated why Chromium atoms 
are the best candidates at the moment for the realisation of such systems. We hope that 
with our work, we have drawn the attention of the statistical mechanics community de- 
voted to the study of long-range interactions to this novel kind of systems and that this 
will be source of inspiration for future collaborations. 
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